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Abstract. A dispersion relation for electromagnetic wave propagation in a strongly 
magnetized cold plasma is deduced, taking photon-photon scattering into account. It 
is shown that the combined plasma and quantum electrodynamic effect is important 
for understanding the mode-structures in magnetar and pulsar atmospheres. The 
implications of our results are discussed. 

PACS numbers: 52.25.Xz (Magnetized plasmas), 52.35.Mw (Nonlinear phenomena), 
52.27.Ep (Electron-positron plasmas) 

1. Introduction 

The quantum electrodynamical (QED) phenomenon of elastic photon-photon scattering, 
due to the interaction of photons with virtual electron-positron pairs, has recently 
received increased attention [1-11]. Several papers are motivated by the desire to 
detect photon-photon scattering in laboratories [1—4], whereas others [5-7] concern 
phenomena that might be relevant when the laser power is further increased to produce 
electric fields strengths close to the Schwinger field ~ 10 18 Vm _1 [7]. Up to now, 
however, observable effects of photon-photon scattering are likely to occur only for 
astrophysical systems [9,10,12,15-18], where the large magnetic field strengths in pulsar 
and magnetar environments [12-14] open up for QED processes to play an important 
role, leading to phenomena such as frequency down-shifting [15, 16] and lensing [17]. 
The frequency down-shifting is a result of so called photon splitting [15,16], which 
is one of the consequences of elastic photon-photon scattering, and the process may 
even be responsible for the radio silence of magnetars [18]. Another QED-process of 
interest in pulsar and magnetar environments is pair-production [19] due to the strong 
field interactions, which lead to the presence of an electron-positron pair plasma in 
the pulsar and magnetar atmospheres. However, with a few exceptions (e.g. [8, 10]), 
we note that most papers considering photon-photon scattering have omitted plasma 
effects when considering electromagnetic wave propagation under these conditions. 

In the present paper we will consider electromagnetic wave propagation at an 
arbitrary angle to a strong external magnetic field B , and include the QED-effects 
associated with that field, as well as the influence of an electron-positron pair plasma. 
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The former effect is described within the framework of the Heisenberg-Euler Lagrangian, 
which constitutes an effective theory of photon-photon scattering [20,21], and the 
latter contribution follows from elementary plasma theory. A comparatively general 
dispersion relation will be derived. It reduces to previous results in a number of 
limiting cases [10,16,22]. In order to determine the contribution from the pair-plasma 
on the propagation properties in pulsar and magnetar atmospheres, we adopt the 
Goldreich- Julian expression for the plasma density [19], and evaluate the dispersion 
relation for field strengths in the pulsar and magnetar range, B ~ 10 8 — 10 10 T. In 
the radio-wave regime it then turns out that for one of the EM-wave polarizations, 
the plasma effects are typically negligible as compared to the QED-effects, whereas 
for the other polarization, the opposite is true in most cases. Noting that important 
processes in pulsar and magnetar environments, e.g. photon splitting, typically involve 
both EM-wave polarizations, we will conclude that QED and plasma effects should be 
simultaneously included when studying radio wave propagation in such environments. 



2. Derivations 



If vacuum fluctuations are taken into account, such as under highly energetic conditions 
(e.g. pulsar plasmas and the next generation of laser-plasma systems), Maxwell's 
equation will be altered by the quantum vacuum self-interaction through the polarization 

P = 2k€ 2 [2(E 2 - c 2 B 2 )E + 7c 2 (E • B)B] (1) 

and magnetization 

M = 2Ke 2 c 2 [-2(E 2 - c 2 B 2 )B + 7(E • B)E] (2) 

respectively, see e.g. [3]. Here k = (a/90ir)(l/e E 2 rit ) gives the strength of the quantum 
vacuum nonlinearity, a ~ 1/137 is the fine-structure constant, E crit = m 2 c 3 /eh ~ 
10 18 Vm" 1 is the Schwinger critical field, m is the electron rest mass, c is the speed of 
light in vacuum, e is the magnitude of the electron charge, and h is Planck's constant 
divided by 2tt. These corrections to Maxwell's vacuum equations are valid as long as 
|E| <C E cr[t and uo <C u e = mc 2 /h, where uo e ~ 8 x 10 20 rads^ 1 is the Compton frequency. 

Next we Fourier decompose the electromagnetic perturbations, which have 
frequencies u and wavevectors k. Maxwell's equations together with the plasma 
equations of motion then yield 

A ab 5E b = 0. (3) 

using index notation. Here the matrix A ab = n a n b — n 2 5 ab + e ab , where n a = k a c/uj, 
n = kcjoj is the plasma refractive index, where k = |k|, e ab = classical + e QED i s the 
dielectric tensor, 

(5) 



e QED = -4£ 
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s denotes the plasma particle species, uj ps = (q 2 n s / 'e m s ) 1//2 is the plasma frequency for 
species s, £ = n 
b a = B^/Bq, and 



species s, £ = k€ c 2 Bq = (a/90Tr)(cB /E crit ) 2 is the dimensionless QED parameter 



1 = -iujd ab + uj cs r ] abl b jl (6) 

with the cyclotron frequency uj cs = q s B /m s for species s, 5 ab is the Kronecker dela and 
rjabc is the totally anti-symmetric unit tensor. From the definition (JHJ) we obtain 



a 
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IUJ 



id 2 — CO 2 
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and the dielectric tensor (jlj) is thus 

,2 



-classical 



ab 



E 



UJ. 



p.s 



vJ A 



(5 ab - b a b b ) + 
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(8) 



We note that the full dielectric tensor depends on the wavevector through the QED 
contribution e^ ED . Freely propagating waves are characterized by the vanishing of the 
dispersion relation D(uj,h) = det(A ab ). Writing the QED-tensor eQ ED in matrix form, 
letting the k-vector lie in the xz-plane, we then obtain 

,2 n \ 



/ 1 



Jib 



-QED 



-4e 





n 2 — 2n 2 ± 




n±n\\ 






\ U ± n \\ ■ 2 

where n\\ = k\\c/u, n± = k±c/u> and the k-vector is written as k = fej_x + fcyz. From 
the classical contributions to A ab is 
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The determinant of the sum of the matrixes (JHJ) and (fTU|) is then evaluated to give the 
dispersion relation 

,2 




x 



UJ 



UJ 



1 + 10e -2J "f - »i (1 - 4 0|- nlnf (1 - 40 



^ uj{u 2 - ujI) 




The dispersion relation (|TT|) is the main result of the present paper. It describes 
wave propagation at any angle to the external magnetic field in a multi-component 
plasma, and it includes the QED effects associated with the external magnetic field. 
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Figure 1. The normalized frequency Q — lu/lo p plotted as a function of the normalized 
wavenumber K = ck/u) p for three cases where k z = 0. The full line represents the 
vacuum dispersion curve, the dashed line represents the strong magnetic held QED 
corrected vacuum dispersion, and the dashed-dotted line depicts the combined effect 
of plasma dispersion and QED effects. The QED corrections corresponds to a magnetic 
held strength of the order of 10 10 T, i.e. held strengths expected around magnetars. 

Thus, it applies to high frequency electromagnetic waves of any polarization, as well as 
electrostatic oscillations and low frequency waves, such as Alfven waves. As a specific 
example of how the plasma dispersion relation is affected by the QED effects we consider 
the case of an electron-positron plasma with uj ~ u p <C |a; c |. The dispersion relation 
relation for the ordinary mode propagating perpendicular to the background magnetic 
field, with strength ~ 10 10 T is depicted in figure 1. A number of limiting cases of 
(JTTJ) have previously appeared in the literature. First, neglecting the QED-effects (i.e. 
letting £ — > 0), we immediately obtain the standard dispersion relation for a cold multi- 
component plasma (see e.g. [22]). Alternatively, letting u p — > 0, we note that the 
dispersion relation depends on the propagation angle relative to the magnetic field. 
Furthermore, we note that the indicies of refraction depend on the polarization even 
without a plasma. These QED-effects due to a strong external magnetic field are 
wellknown (often referred to as "birefringence of vacuum"). Our dispersion relation 
in the limit u p — > agrees with those of previous works, see e.g. [15,16]. The combined 
contribution from the QED-effects due to a strong magnetic field and a non-zero plasma 
density have previously been considered [10] in the limit of parallel propagation and 
allowing for large amplitudes. Taking the limit of a small wave amplitudes in the 
dispersion relation (11) of reference [10], and letting n± — > in (jll)) we obtain agreement 
with [10]. 
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3. Conclusion 

QED-effects associated with the external magnetic field are likely to be of importance in 
environments with extreme magnetic fields, in particular in the vicinity of astrophysical 
objects like pulsars and magnetars. For example, the radio silence of magnetars is 
assumed to be connected with QED-effects associated with the magnetar fields [18], 
which could reach 10 10 — 10 11 T close to the surface. However, in the same environments, 
we also expect the presence of an electron-positron plasma [19]. Thus we evaluate 
(fTTj) with J2 S = Sep' w here e and p denotes electrons and positrons, respectively. 
Considering propagation at an arbitrary angle to the magnetic field in an electron- 
positron plasma, letting uj P£jPP ~ u> <C |^ce,cp| , using ( < 1 and noting that the 
factor [^2 ep uj ps uj cs /uj(uj 2 — u 2 s )] 2 then becomes negligibly small (due to the approximate 
cancellation of the electron and positron contributions), we find from (fTTj) that the 
dispersion relation separates in two modes that can be approximated by 

,2 



LOZ. 



1 -" 2+8 ^ + ;^r^r (12) 



and 



oj 2 



(1 - n 2 )(l - 40 - (-14£ + J ) (1 - nf) « 0, (13) 

where tu p = (to pe + uj 2 ^ 1 ^ 2 is the total plasma frequency, and tu c = eB /m is the 
magnitude of the electron (or positron) cyclotron frequency. In the vicinity of pulsars 
or magnetars where u c ~ 10 19 — 10 21 rads^ 1 , the last term of (|12j) is negligible unless the 
plasma density is extremely high. Omitting that term, the dispersion relation (JT2")l is 
then the same as that used in [16] for the high phase velocity mode when considering 
photon-splitting. Similarly the ordinary mode described by (|13jh reduces to the mode 
with the lower phase velocity of reference [16] when the plasma is removed. However, for 
the latter dispersion relation we note that a relatively modest plasma density is enough 
to significantly affect the propagation properties in the radio wave regime. To make a 
concrete estimate, we adopt the Goldreich-Julian density 

ll5 / 0.1 \ / -Bp U i sar \ _g 



, iGJ = 7 x 10" m- (14) 

where r is the pulsar period time (in seconds) and -B pu i S ar the pulsar magnetic field (in 
tesla). The pair plasma density is expected to satisfy n e = n p = Mn G j, where M is 
the multiplicity [19,23]. Moderate estimates then give M — 10 [23]. Choosing this 
value and letting r = 1 s, we note that for magnetar field strengths, -B pu isar = 10 10 T, 
the term due to the plasma oc uj 2 /u 2 in (fTB*)) dominates over the term due to QED 
oc 14£ for frequencies up to uj ~ 10 14 — 10 15 rads _1 , i.e. in the infrared regime and 
below. Furthermore, we note that photon splitting [15, 16] as described by standard 
QED (i.e. with zero plasma density) requires that the phase velocity of the dispersion 
relation in (|12jl is higher than that of 1)13)1 . While this is always true in the absence 
of a plasma, we note that for wave frequencies in the infra-red regime and below, the 
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Goldreich- Julian density given by (J 14)) is enough to increase the phase velocity of the 
mode in (JT^) above that of (|T2|) . unless we choose the period time r extremely low. 
Thus we conclude that photon-photon splitting as described by vacuum theories is not 
likely to apply to magnetar atmospheres, unless the pair-production [19] responsible for 
the Goldreich- Julian expression is effectively suppressed. Wave cascade processes as a 
mechanism to explain the radio silence of magnetars [18] could still be possible, but for 
densities of the order of ([T4"|) . plasma nonlinearities are likely to dominate over the pure 
QED effects. 
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